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Abstract — This paper presents development of prototype prod- 
ucts for terrestrial ecosystems in preparation for the future imag- 
ing spectrometer planned for the Hyperspectral Infrared Imager 
(HyspIRI) mission. We present a successful demonstration ex- 
ample in a coniferous forest of two product prototypes: fraction 
of photosynthetically active radiation (PAR) absorbed by chloro- 
phyll of a canopy (fAPAR ch i) and leaf water content (LWC), 
for future HyspIRI implementation at 60-m spatial resolution. 
For this, we used existing 30-m resolution imaging spectrometer 
data available from the Earth Observing One (EO-1) Hyperion 
satellite to simulate and prototype the level one radiometrically 
corrected radiance (L1R) images expected from the HyspIRI 
visible through shortwave infrared spectrometer. The HyspIRI- 
like images were atmospherically corrected to obtain surface re- 
flectance and spectrally resampled to produce 60-m reflectance 
images for wavelength regions that were comparable to all seven 
of the MODerate resolution Imaging Spectroradiometer (MODIS) 
land bands. Thus, we developed MODIS -like surface reflectance in 
seven spectral bands at the HyspIRI-like spatial scale, which was 
utilized to derive fAPAR ch i and LWC with a coupled canopy-leaf 
radiative transfer model (PROSAIL2) for the coniferous forest. 
With this paper, we provide additional evidence that the fAPAR ch i 
product is more realistic in describing the physiologically active 
canopy than the traditional fAPAR parameter for the whole 
canopy (fAPAR canopy ), and thus, it should replace it in ecosystem 
process models to reduce uncertainties in terrestrial carbon cycle 
and ecosystem studies. 
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Index Terms — Earth Observing One (EO-1) Hyperion, 
fAPAR can0 p y , fAPAR c hi, foliar moisture content, Hyperspectral 
Infrared Imager (HyspIRI), leaf water content (LWC), terrestrial 
carbon cycle. 

I. Introduction 

T HE HYPERSPECTRAL Infrared Imager (HyspIRI) 
mission was described by the National Research Coun- 
cil in its Decadal Survey Report (http://www.nap.edu/ 
catalog.php?record_id=11820) to address terrestrial ecosystem 
science as one of the next-generation NASA satellite missions. 
The HyspIRI mission is envisioned to carry two spectral in- 
struments, both with ground spatial resolutions of 60 m — a 
visible to shortwave infrared (VSWIR) continuous spectrum 
hyperspectral imager (10-nm spectral sampling) and a multi- 
channel thermal infrared (TIR) imager. The Earth Observing 
One (EO-1) Hyperion (launched in November 2000) is still 
operating and serves as the heritage satellite instrument for 
HyspIRI’s VSWIR spectrometer, but it only captures 7.5-km- 
wide ground strips, and its 30-m resolution images are acquired 
through user/system requests. In contrast, HyspIRI will be a 
global survey mission, and its VSWIR instrument will have 
60-m pixels across a 150-km-wide ground swath, collected on 
an equatorial 19-day repeat cycle. Consequently, existing Hy- 
perion data provide an excellent tool for product development 
in anticipation of the HyspIRI and other spacebome imaging 
spectrometer missions. 

The absorbed photosynthetically active radiation (APAR) 
fraction for a whole vegetation canopy (fAPAR can0 py, also 
denoted as FAPAR or FPAR [1]— [3]; see Appendix A for 
equations) is an essential climate variable [4]— [7] needed 
to estimate and monitor vegetation productivity on a global 
basis. However, fAPAR canopy includes both photosynthetic 
and nonphotosynthetic components and has not provided 
consistent relationships to photo synthetic processes at the 
ecosystem scale [8]— [1 1]. This is because the APAR available 
to support vegetation photosynthesis (APARpsn) is typically 
overestimated by fAPAR canopy . However, the APAR fraction 
associated with the chlorophyll-containing component 
(fAPAR c hi, (A3) in Appendix A) consistently and correctly 
represents the physiologically active photosynthetic sector of 
the canopy under optimal (e.g., fully green) and less optimal 
(e.g., mixtures of green and senescent vegetation) conditions 
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affecting physiological responses. In other words, 
fAPARp S N = fAPAR ch i [11]. 

We recently demonstrated that fAPAR c hi is superior to the 
use of fAPAR canopy in model simulations with gross primary 
production (GPP) or gross ecosystem production to estimate 
light use efficiency (LUE), defined as GPP/APARpsn [11]. The 
fAPAR c hi retrievals were estimated from space for a deciduous 
aspen forest using five of the seven MODerate Resolution 
Imaging Spectroradiometer (MODIS) spectral land bands from 
Collection 4 daily products, for which product quality was 
insufficient for our model retrievals in two land bands (B3, blue; 
B7, SWIR2). Those earlier results were indirectly validated 
by comparing LUE measured in situ at the tower (LUE t0 wer) 
to the LUE determined from our remote sensing/modeling 
approach for the forest’s chlorophyll component (LUE c hi = 
GPP/APAR c hi, where APAR c hi = fAPAR c hi * PAR). LUE c hi 
matched well with LUE t0 wer, while the widely used LUE 
describing the whole canopy (LUE canopy = GPP/APAR canopy , 
where APAR canopy = fAPAR canopy * PAR) did not. Therefore, 
we recommended that fAPAR c hi should replace fAPAR canopy 
to estimate canopy parameters related to photosynthesis for cli- 
mate models and land-atmosphere interaction models [4], [6]. 
However, further evidence should be pursued. 

The spectral range for both the EO-1 Hyperion and future 
HyspIRI VSWIR imaging spectrometers is between 0.4 and 
2.5 /im, which spans the spectral range of the MODIS 1-7 
land bands. In the present study, the fAPAR c hi algorithm that 
was previously developed to ingest five MODIS bands was 
modified to utilize all seven MODIS land bands from the more 
radiometrically rigorous Collection 5 products. We wish to 
know how inclusion of these additional bands and higher spatial 
resolution satellite observations (60 versus 500 m) affected and 
improved retrievals of fAPAR c hi from HyspIRI-like VSWIR 
radiance images simulated from EO-1 Hyperion images. We 
also simultaneously retrieved leaf water content [LWC, (A4)]. 
The Normalized Difference Vegetation Index [NDVI, (A5)] and 
Enhanced Vegetation Index [EVI, (A6)] were used to estimate 
fAPAR canopy and fAPAR c hi, respectively [12]— [15], and the 
Land Surface Water Index [LSWI, (A7)] was used to estimate 
foliar moisture content [14]. Note that fAPAR canopy is a linear 
function of NDVI [12]. 

The goal of this paper was to apply the modified fAPAR c hi 
and LWC algorithm to a coniferous forest in a heteroge- 
neous landscape, to demonstrate the advantages of the revised 
algorithm to observations with spectral bands spanning the 
full optical spectrum at much high spatial resolution than is 
possible with MODIS. Our specific objectives were to test 
the following hypotheses: 1) fAPAR c hi and LWC provide 
unique information, as compared to existing indices such as 
EVI, fAPAR can o Py (linear function of NDVI), and LSWI, and 
2) fAPARchi and LWC retrievals benefit from higher spa- 
tial resolution and additional spectral band inputs. We begin 
by describing the approach to obtaining prototype HyspIRI 
VSWIR radiance images and then describe the modification 
of the fAPARchi -LWC algorithm from five- to seven-band ver- 
sions. Next, we present the HyspIRI outputs of the seven-band 
algorithm and comparisons of our model retrievals for 500- 
versus 60-m pixels, and comparison of 500-m retrievals using 
five versus seven spectral bands, followed by the summary 
conclusion. 


II. Methods 

A. Satellite Image Preprocessing 

1) Spatially Scaling Up the EO-1 Hyperion Radiance Im- 
ages to 60 m: The EO-1 Hyperion images have a spatial 
resolution of 30 m. We spatially scaled up the Hyperion Level 
One radiometrically corrected Radiance (L1R) data to 60 m 
by averaging Hyperion 30-m pixels in four-pixel blocks [16], 
[17] to obtain a spatially relevant prototype of 60-m HyspIRI 
L1R data, also achieving an average signal to noise response 
comparable to that expected for HyspIRI (> 400: 1). These 
measured radiances were divided by solar irradiances above 
the atmosphere to obtain the apparent top-of-atmosphere (TOA) 
reflectances. 

2) Atmospheric Correction With the ATREM Routine: In 
order to use spectral imaging data for quantitative remote 
sensing of land surfaces, the absorption and scattering effects 
of atmospheric gases and aerosols must be removed [18]. The 
HyspIRI-like L1R images (at 60 m) were atmospherically cor- 
rected using an updated version of the ATmosphere REMoval 
Algorithm (ATREM) with which a line-by-line model was 
used to calculate atmospheric gaseous transmittances [19], [20]. 
The surface reflectances were derived from the apparent TOA 
reflectances using the simulated atmospheric gaseous transmit- 
tances and the simulated molecular and aerosol scattering data. 
During retrievals, the integrated water vapor amount on a pixel- 
by -pixel basis can be directly derived from the 0.94- and the 
1.14-/im atmospheric water vapor absorption features, a spe- 
cial advantage conveyed by continuous spectrometer data. The 
transmission spectrum of water vapor (H 2 0), carbon dioxide 
(CO2), ozone (O3), nitrous oxide (N 2 0), carbon monoxide 
(CO), methane (CH 4 ), and oxygen (0 2 ) in the 0.4-2.5-/im 
region was simulated based on the derived water vapor value 
and the solar and observational geometry, and through the use 
of narrowband spectral models. However, the scattering effect 
due to atmospheric molecules and aerosols was determined 
with the 6S computer code [21]. 

3) Spectrally Combining HyspIRI-Like Surface Reflectance 
Bands to Simulate the MODIS Bands 1-7: After obtaining at- 
mospherically corrected HyspIRI-like 60-m surface reflectance 
images, the spectral values were averaged across 3-6 contigu- 
ous 10-nm Hyperion bands within the defined MODIS band 
ranges (Table I) using spectral response functions to obtain 
reflectances spectrally comparable to those from MODIS [22]. 
This yielded an image that was spatially HyspIRI-like but spec- 
trally MODIS -like. To summarize, spatially scaling-up from 
30 to 60 m was performed on the Hyperion L1R radiance 
image; after which, the ATREM atmospheric correction was 
performed on the 60-m HyspIRI-like L1R radiance image (60 x 
60 HyspIRI pixel block) to obtain a 60-m HyspIRI-like surface 
reflectance image, followed by spectral averaging to obtain 
surface reflectance in each of the seven MODIS -like bands 
(Table I). 

The vegetation canopy parameters derived from the 
fAPAR c hi algorithm [11] included the following: leaf internal 
structure (N), leaf dry matter (C m ), leaf water thickness (C w ), 
and leaf pigment content (C a b). In preliminary model runs, 
all seven (of 36) MODIS land bands (1-7) were found to be 
sensitive to N and C m , whereas bands 1,3, and 4 were sensitive 
to C a b, and bands 5, 6, and 7 were sensitive to C w [23]. 
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TABLE I 

Spectral Ranges Covered by the MODIS 
AND HYPERION/HYSPIRI BANDS 


Spectral range/ Band Width 

MODIS band # 

Hyperion/HyspIRI band # 

459 - 479 nm/ 20 nm 

3 (blue) 

11-13 

545 - 565 nm/ 20 nm 

4 (green) 

20-22 

620 - 670 nm/ 50 nm 

1 (red) 

27-32 

841 -875 nm/ 34 nm 

2 (NIRO 

49-52 

1230- 1250 nm/ 20 nm 

5 (NIR 2 ) 

108-111 

1628 - 1652 nm/ 24 nm 

6 (SWIR0 

148-150 

2105- 2155 nm/ 50 nm 

7 (SWIR0 

195-200 


B. Algorithm to Derive fAPAR chl and LWC Using PROSAIL2 
and the Metropolis Approach 

A complete description of the PROSAIL2 model and 
Metropolis approach, as applied to five spectral MODIS land 
bands, is given in a recent publication [11]. The coupled 
canopy-leaf radiative transfer model utilized in this PROSAIL2 
algorithm is based on the SAIL2 canopy radiative transfer 
model and the PROSPECT leaf radiative transfer model. Here, 
we provide an overview and highlight the changes introduced 
in the revised approach. Additional details and information are 
provided in Appendix A. 

In brief, a vegetation canopy can be partitioned into leaf and 
nonleaf (referred to as stem) components. A leaf can be fur- 
ther partitioned into chlorophyll, nonphoto synthetic pigments 
(referred to as brown pigment Cbrown), water, and dry matter. 
The PROSAIL2 model has 14 biophysical and biochemical 
variables (see Appendix A), including five leaf variables that 
simulate leaf optical properties (N, C a b, C m , C w , and CbrownX 
a soil/litter variable that simulates soil/litter optical properties 
(SOILa), and a variable that simulates stem optical properties 
(STEMa). 

We modified the previous MODIS fAPAR c hi algorithm [11] 
by replacing the five-band likelihood function with that for 
seven-band surface reflectances obtained from Section II-A3 
[(1) and (2)]. The Markov Chain Monte Carlo (MCMC) method 
(Metropolis) is employed for inversion. This method assumes 
that the observed spectral reflectances 0Q = [xn , . . . ,Xi P ] f 
(p = 7; Table I) differ from the model predicted values Ui = 
[uii, . . . , Ui p \' according to a mean zero p-variate Gaussian 
error model that results in the likelihood function 




1 c -(X i -U i )'Y,- 1 (X i -U i )/2 

(v^7r) p |£|! 


(l) 


where n is the number of data points sampled according to 
equation (8) of Zhang et al. [11] and E is the 
variance-covariance matrix of X 


_ ( s ij)pxp 
e 

1 U 

Sij = ^ ^ ( %ki ^ ki ) {p^kj ^ kj ) •> ^ 5 j = 1 •> • • • 5 P- (2) 

k=l 


The new seven-band fAPAR c hi algorithm provides simulta- 
neous solutions for C w and C m , enabling the solution of LWC 
[see (A4)]. The full solution of the 14 parameters is a statistical 
posterior distribution based on the radiative transfer model and 
the remote sensing observation (see Section II- A3). 

The 14 parameters and the derived fAPAR c hi and LWC may 
be grouped into three classes based on their posterior statistical 
distributions: well-constrained, edge-hitting, and poorly con- 
strained. The posterior statistical distributions can provide the 
mode(s) of the 14 variables, fAPAR ch i, and LWC (if they exist) 
(please see Zhang et al. [24] for details). A mode for a variable 
is one “traditional” (best) point solution, i.e., the most likely 
value of the variable to fit both the PROSAIL2 model and the 
remote sensing observation. Lrom the case study of this paper 
(see Section II-C for site description), we discovered that there 
was one and only one mode for fAPAR c hi and for LWC per 
satellite observation. The reason is that chlorophyll, leaf water, 
and dry matter have unique spectral characteristics. These 
components can be distinguished using spectral information of 
the seven bands, and they will not mess up with each other or 
other components of a canopy. 

The EVI [25], NDVI [26], fAPAR canopy , and LSWI [14] 
were also calculated (see Appendix for equations) and com- 
pared with results for fAPAR c hi or LWC, as appropriate, at both 
MODIS (500 m) and HyspIRI (60 m) spatial resolutions. We 
also calculated fAPAR canopy for the whole canopy based on 
a widely used formula which relates fAPAR canopy and NDVI 
(e.g., [12], [13], and [15]) 

fAPAR canopy = 1.24 x NDVI - 0.168. (3) 

We also utilized the same seven-band approach [(1) and (2)] 
to retrieve fAPAR c hi and LWC from MOD09A1 [the eight-day 
composite reflectance MODIS product (M)] acquired on day 
185 (July 3, 2008), which was close to the acquisition date of 
the original Hyperion image (June 28, 2008). This enabled us 
to compare the 60-m HyspIRI-like product (H) with MODIS 
500-m product. 


C. Study Site 

The study site (Lig. 1) was a Douglas fir forest surrounding 
an instrumented tower (hereafter DL49: 49° 52' N, 125° 20' W, 
300-m elevation) in the Canadian Carbon Program network. 
The DL49 is located at the eastern side of Vancouver Island, 
BC, Canada, and the forest stand around the tower (indicated 
by a circle in Pig. 1) is mainly comprised of Douglas fir with 
some western red cedar and western hemlock [27]. The study 
area was a 120 x 120 grid formed by Hyperion pixels at 
the original 30-m spatial resolution around the DP49 site. A 
true color red/green/blue (RGB) image at the 60-m HyspIRI 
pixel spatial resolution (Pig. 1) can be compared with the land 
cover map produced using the ISODATA method of ENVI 
(Pig. 2), which utilized the surface reflectances of the original 
HyspIRI bands (Table I) for land cover classification. “Unveg- 
etated” areas are associated with roads or sparse vegetation. 
Harvested areas show various stages of forest regeneration. 
Wetter forested areas are dominated by hemlock, alder, and 
maple (personal communication, Nicholas Coops, University of 
British Columbia, Vancouver, BC, Canada). 
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Fig. 1. (a) Location of the Douglas fir site (DF49) on Vancouver Island, BC, Canada, (b) True color RGB image for the DF49 area using simulated HyspIRI data 

on DOY 180, 2008 (June 28, 2008), where the circle designates the fetch of the DF49 flux tower. 



Legend 


Unvegeteted area 
Harvested area 
Douglas-fir 
Wetter forest area 


Fig. 2. Land cover map for the DF49 area using simulated HyspIRI data based 
on the EO-1 Hyperion image collected on DOY 180, 2008 (June 28, 2008). 


III. Results 
A. HyspIRI-Like Results 

Here, we present the results of the revised seven-band 
fAPAR c hi-LWC algorithm (described in Section II-B) applied 
to the mid-summer HyspIRI-like L1R radiance image in the 
vicinity of the DF49 Douglas fir tower site (described in 
Section II-C) for the purpose of developing and evaluating 
prototype products. 

Fig. 3 has three subfigures for our study site showing 
the spatial distributions for three of the variables of interest: 
1) fAPAR c hi; 2) fAPAR canopy [based on NDVI using (3)]; 
and 3) EVI, all presented on the same relative scale between 
0.0 and 1.0. Clearly, the values in the fAPAR c hi map are 
substantially and statistically lower than those exhibited by the 
fAPAR canopy map. Table II lists the mode, mean, and median 
values computed for the study area shown in Figs. 1(b), 2, 
and 3(a)-(c) for fAPAR ch i, EVI, fAPAR canopy , and NDVI. 
These statistical values reveal that fAPAR c hi and EVI provide 
substantially lower values than NDVI and fAPAR canopy . 

The vegetation photosynthesis model (VPM) [14] assumes 
fAPAR c hi = EVI, but we find that this assumption is not always 
correct. Moreover, overestimates of APARpsn result if we 
assume fAPAR c hi = fAPAR canopy . 

The LWC and LSWI maps are shown in Fig. 4. The LWC 
values of wetter forest areas (0.600-0.909) differ substantially 


from LSWI values (0.202-0.483). Although the LSWI has been 
shown to represent water status of vegetation in some studies 
[14], [28], this index cannot distinguish canopy water from 
background water (e.g., soil water). The LWC of different 
plant species among land cover types might vary, as shown 
in Fig. 4(a) for the LWC dynamics per class. The broad-leaf 
deciduous leaves in the wetter forest areas had higher average 
LWC than the Douglas fir leaves [Fig. 4(a)]. 

The histograms for fAPAR c hi, fAPAR canopy , and EVI and 
the histograms for LWC and LSWI are shown in Figs. 5 and 6, 
respectively. Peak frequencies occurred at very different values: 
fAPAR ch i (0.559), EVI (0.437), and fAPAR ca nopy (0.893), but 
peak frequencies for LWC and LSWI occurred at similar value 
(^0.49). The fAPAR c hi parameter displayed minor mode at 
0.714. For forests, the magnitude is higher, and the range is 
wider for fAPAR c hi as compared to EVI [Fig. 3(a) and (c)]. 
The frequency of the mode derived from the HyspIRI-like LWC 
map is about three times that of the comparable mode for the 
LSWI (Fig. 6). Scatter plot comparisons are also shown for 
these pairs: fAPAR c hi versus fAPAR canopy , fAPAR c hi versus 
EVI, and LWC versus LSWI (Fig. 7). The first pair exhibits 
that fAPARcanopy is greater than fAPAR c hi. While values are 
closer between EVI and fAPAR c hi, the slope relating these 
two variables clearly deviates from the 1:1 line. No apparent 
correlation is seen for the third pair (LWC and LSWI). LSWI 
cannot be used to replace or predict LWC well, as shown in 
Fig. 7(c). 

The pixels classified as “unvegetated” were recently har- 
vested (personal communication, Nicholas Coops, University 
of British Columbia). The fAPAR c hi values for those pixels are 
close to zero, as should be expected for areas without green 
vegetation. That is to say that fAPAR c hi has a physical and 
physiological meaning. However, the NDVI and EVI values of 
those pixels, which are greater than 0.4 and 0.2, respectively, 
indicate the presence of some green vegetation. NDVI (and the 
derived parameter fAPAR can opy) saturates for pixels with a leaf 
area index greater than three [25], while fAPAR c hi does not. 


B. Comparing HyspIRI-Like and MODIS Parameter Estimates 

From the processed MODIS image (July 3, 2008; DOY 185), 
we selected the pixel that covers the DF49 site. The HyspIRI- 
like maps have considerably more spatial details than those 
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Fig. 3. HyspIRI-like maps for the DF49 area, (a) fAPAR c hi. (b) fAPAR can opy computed from NDVI (3). (c) EVI. Data were simulated from the mid-summer 
Hyperion image acquired on DOY 180, 2008 (June 28, 2008). 


TABLE II 

Values of Mode(s), Mean, and Median of the fAPAR chh EVI, 
fAPAR canop y, AND NDVI VALUES FOR THE AREA SHOWN 
in Fig. 3. Number of HyspIRI-Like Pixels = 3600 



mode(s) 

mean 

median 

fAPAR c hi 

0.559, 0.714 

0.520 

0.544 

EVI 

0.437 

0.460 

0.452 

fAP AR can0 py 

0.893 

0.824 

0.870 

NDVI ** 

0.855 

0.800 

0.837 


map not shown. 


based on the MODIS image. To compare our retrievals from 
HyspIRI-like images with those based on the MODIS image, 
using the same seven-band approach [(1) and (2)], we chose 
only the HyspIRI pixels that fall within 240 m of the DF49 
tower site. Forty-nine HyspIRI pixels were selected (Table III). 

In the comparison of the seven-band HyspIRI-like (H) and 
true MODIS (M) products (Table III), the mean fAPAR ch i val- 
ues (H) and the single MODIS fAPAR c hi value (M) were fairly 
similar. The mean H and M satellite values for fAPAR canopy 
were comparable to the tower-based canopy-level fAPAR (0.94 
[29]) determined from the DF49 tower radiation measurements 
for the same period. The green leaf fAPAR (0.79 [30]) esti- 
mated using the approach developed by Chen et al. [31], [32] 
represented the combined effects of chlorophyll fAPAR and the 
fAPAR of leaf dry matter and brown pigments of the canopy for 
the same period. Thus, the green leaf fAPAR was intermediate 
between estimates for fAPAR c hi and fAPAR can0 py. 


LWC provides quantitative information on foliar moisture 
content. It is not only a critical indicator of vegetation growth 
status but also an important factor in the canopy susceptibility 
to the fire ignition process. Our LWC retrievals for H and M 
(Table III) both fall within the published Douglas fir LWC range 
(between 0.67 and 0.44), which includes both young and old 
leaves [33], [34]. The mean LWC value for the 49 HyspIRI-like 
pixels [Fig. 4(a)] is the same as that for the MODIS single pixel 
LWC (-0.49; Table III). 

In addition to application of the seven-band algorithm to the 
MODIS pixel that covers the DF49 site, we also calculated 
fAPAR c hi and LWC using the original five-band algorithm [11]. 
The seven-band versus five-band fAPAR c hi histograms of the 
pixel have the same mode value (0.533), with only slightly 
different standard deviations (0.071 versus 0.072). However, 
while the seven-band versus five-band LWC histograms at the 
tower site have the same mode value, their standard deviations 
differ (0. 169 versus 0. 174), such that less uncertainty is incurred 
using the seven-band version. 

IV. Discussion 

This paper describes how to estimate two products 
(fAPAR c hi and LWC) for HyspIRI L1R radiance and presents 
some initial prototype results. In addition to our primary prod- 
ucts (fAPARchi and LWC), we also determined values and 
map products for the EVI, fAPAR can opy, and LSWI. Although 
fAPARchi and EVI were the most similar, the range of values 
and the modes for fAPAR c hi were larger than those obtained 
for the EVI (Fig. 5(a) and (c); Table II). Likewise, the dy- 
namic range for fAPAR can opy was smaller than the range 
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(a) LWC (b) LSWI 



Fig. 4. HyspIRI-like maps for the DF49 area, (a) LWC. (b) LSWI. Data were simulated from the mid-summer Hyperion image acquired on DOY 180, 2008 
(June 28, 2008). 



Fig. 5. Histograms for (a) fAPAR c hi, (b) fAPAR can0 py, and (c) EVI for the DF49 area shown in the Fig. 3 maps. 



(b) 


025 
0.2 
0 15 
0 1 
0.05 


0 "0 2 0.4 0.6 0,0 


Fig. 6. Histograms for (a) LWC and (b) LSWI for the DF49 area shown in the Fig. 4 maps. 


i 



fAPAR^i 


fAPARchi 


LWC 


Fig. 7. Comparisons for (a) fAPAR c hi versus f APAR can opy , (b) fAPAR ch i versus EVI, and (c) LWC versus LSWI. Data are derived from the simulated HyspIRI 
image acquired on DOY 180, 2008 (June 28, 2008). 
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TABLE III 

Comparison at the DF49 Tower Site of the fAPAR ch i, 
f APAR canop y , NDVI, EVI, LSWI, and LWC Values From the 
Simulated HyspIRI-Like (60 m) Image With the MODIS Image 
(500 m) and Published Field Measurements (H: HyspIRI-Like 
Data; M: Seven-Band MODIS-Based Data) 



Parameters 

Value 
(± STDEV) 

Sample size, 
Comments 

fAPAR 

related 

parameters 

fAPAR C hi 

H 

0.583 ± 0.038 

n = 49 pixels 

M 

0.533 

single pixel 
value 

fAP AR C anopy 

H 

0.907 ±0.006 

n = 49 pixels 

M 

0.949 

single pixel 
value 

Tower based 

0.94 

Hember et al. 
[29] 

day 180-185 
of 2008 

NDVI 

H 

0.867 ±0.005 

n = 49 pixels 

M 

0.901 

single pixel 
value 

EVI 

H 

0.449 ±0.017 

n = 49 pixels 

M 

0.445 

single pixel 
value 

fAPARg reen _| ea f 

Tower based 

0.79 

Hilker et al. 
[30] 

day 180-185 
of 2008 

Leaf 

water 

related 

parameters 

LWC 

H 

0.494 ± 0.008 

n = 49 pixels 

M 

0.493 

single pixel 
value 

Field 

measurements 

0.44-0.67 

Agee et al. 
[33], 

Keyes [34] 

growing 
season of 
multiyears 

LSWI 

H 

0.506 ±0.013 

n = 49 pixels 

M 

0.510 

single pixel 
value 


for fAPAR c hi values [Fig. 5(a) and (b)]. When comparing 
fAPAR c hi and EVI, we found that, when fAPAR c hi = 0.5, the 
EVI range is 0.384-0.533, but when EVI = 0.5, the fAPAR ch i 
range is 0.410-0.686 [Fig. 7(b)]. 

In addition to better spatial detail, one advantage of the 
seven-band fAPAR c hi and LWC algorithm is that it does not 
need land cover type information as an input to run the 
model inversion, whereas the MODIS standard fAPAR ca nopy 
(i.e., FPAR) product does. The seven-band algorithm can pro- 
vide fAPARchi and LWC products with less uncertainty (e.g., 
smaller standard deviations) as compared to results obtained 
with the previous five-band algorithm, even for a relatively 
homogeneous forest area [circle in Fig. 1(b)]. The outputs of the 
algorithm (fAPAR c hi and LWC) can be used for seasonal anal- 


ysis, interannual analysis, phenological study, and land use and 
land cover change research — including disturbance studies and 
disaster monitoring (e.g., fire, drought, and flooding). Our study 
also demonstrates the flexibility that an imaging spectrometer 
allows for interinstrument comparisons. 

The remote sensing community uses three groups of in- 
version strategies: MCMC approaches, look-up tables, and 
gradient-based approaches. With the Metropolis approach (an 
MCMC method), we can globally search for the optimal solu- 
tion, a posterior distribution. However, look-up-table methods 
provide fixed step lengths for all parameters before inversion, 
whereas gradient-based methods can only search local optima 
and rely on initial guesses. We anticipate that application of 
our algorithm to satellite images will be useful for the current 
and future national and international research projects that rely 
on remotely sensed data, including the North American Carbon 
Program. 


V. Conclusion 

We have successfully demonstrated here that the two prod- 
ucts (fAPAR c hi and LWC) provide unique information. The 
most important finding in this paper is that fAPAR c hi values 
differ from those for EVI, NDVI, and fAPAR canopy in most 
cases. EVI is not always equal to fAPAR c hi. We have also 
found the following: fAPAR ch i 7 ^ fAPAR canopy (or NDVI), and 
LWC 7 ^ LSWI. In other words, we have rejected the null hy- 
potheses that equate the EVI and fAPAR canopy with fAPAR c hi 
or LSWI with LWC. HyspIRI also has the potential to provide 
the spatial variance of the two products that cannot be extracted 
from MODIS. We have realized that real HyspIRI images, or 
those of another future imaging spectrometer, will differ in 
some ways from Hyperion, which can be taken into account 
with the at-launch version of the algorithm. 


Appendix 

In brief, a vegetation canopy can be partitioned into leaf 
and nonleaf (hereafter referred to as stem) components. A leaf 
can be further partitioned into chlorophyll, nonphoto synthetic 
pigments (hereafter referred to as brown pigment Cb row n), 
water, and dry matter (or C a b, Cbrown, C w , and C m ). The 
PROSAIL2 model has 14 biophysical and biochemical vari- 
ables: plant area index, stem fraction (SFRAC), cover fraction, 
stem inclination angle (STINC), stem BRDF effect variable 
(STHOT), leaf inclination angle (LFINC), leaf BRDF effect 
variable (LFHOT), five leaf variables that simulate leaf optical 
properties (N, C a b, C m , C w , and CbrownX one soil/litter variable 
that simulates soil/litter optical properties (SOILa), and one 
variable that simulates stem optical properties (STEMa). The 
MCMC method (Methopolis) is employed for inversion. 

One can calculate fAPAR canopy [ 12 ] and fAPAR c hi [11] with 


APAR canopy = APAR ch \ + APAR dry 

matter 


fAPAR, 


AP AR hrmvn pigment 
+ AP A._R s tem 

APAR C 


‘'Canopy 


canopy 


PAR 0 


fAPAR C h\ = 


APAR C h\ 

PARa 


(Al) 

(A2) 

(A3) 
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where PARq is the incoming PAR at the top of the canopy and 
APAR can0 py, APARchl, APARdry matter? APARforown pigment? 
and APARgtem are PAR absorbed by canopy, chlorophyll in 
leaf, dry matter in leaf, brown pigment in leaf, and stem, 
respectively. One has to know the value of PAR 0 and the values 
of the 14 parameters to calculate APAR c hi and APAR can0 py 
in (Al). One may assume PARo be any positive value to 
calculate fAPAR c hi and fAPAR canop y because they are ratios 
[(A2) and (A3)]. We present NDVI-based fAPAR can0 py for this 
paper because of the linear relationship between fAPAR canop y 
and NDVI (3), which is based on the simulation study using the 
SAIL [12]. 

The leaf water thickness (C w in grams per square centimeter 
or centimeter) and leaf dry matter (C m in grams per square 
centimeter) are two of the 14 parameters of PROSAIL2, and 
the inversion algorithm provides their posterior distributions as 
outputs. LWC is defined as the fraction of leaf water weight to 
fresh leaf weight [35]. That is to say 


LWC = —^~ 

U VU + U? 


(A4) 


EVI [25], NDVI [26], and LSWI [14] are also calculated: 


EVI = 2.5 x 


PNIRi —Pre d 

PNIRi + 6.0 X p red - 7.5 X /9 b lu e + 1.0 


N DVI = ^ mRl /9red 

PNIRi +Pred 

LSWI = ^ NIR i PSWIRi 
PSWIRi + PSWIRi 


(A5) 

(A6) 

(A7) 


where p is the reflectance. 
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